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The deconfinement phase transition is studied in the FL model at finite temperature and chemical 
potential. At MFT approximation, the phase transition can only be the first order in the whole fj,— T 
phase plane. By a Landau expansion we further study the phase transition order and the possible 
phase diagram of deconfinement. We discuss the possibilities of second order phase transitions in 
FL model. By our analysis the cubic term in the Landau expansion could be cancelled by the high 
order fluctuations. By an ansatz of the Landau parameters, we obtain the possible phase diagram 
with both first and second order phase transition including the tricritical point which is similar to 
that of the chiral phase transition. 
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^ . I. INTRODUCTION 

' It is generally believed that at sufficiently high temperatures and densities there is a QCD phase transition from 
Or normal nuclear matter to QGP p], [|J ■ Theoretically there are two kinds of phase transitions associated with different 
O ■ symmetries for two opposite quark mass limit. For Nf = 2 + 1 massless quark flavors, the QCD lagrangian posses 
a chiral symmetry SU(Nf)n x SU(Nf)L, which is associated with the chiral phase transition. In the heavy quark 
. limit, QCD reduces to a pure SU(N C ) gauge theory which is invariant under a global Z(N C ) center symmetry. This 
symmetry is associated with the deconfinement phase transition. The orders of these phase transitions have been 
y—i , studied extensively [H-IH and still remained to be an interesting problem For chiral phase transition at finite 

' temperature in the chiral limit, the quark-antiquark condensate (qn.QL) serves as a good order parameter. The order 
of the phase transition depends on the quark flavors. For Nf = 3 massless quark flavors, it is a first order phase 
transition. For Nf — 2 massless quark flavors, it is a second order phase transition. At finite densities, the chiral phase 
transition have been studied by many effective models (ToMT^ j . It is generally regarded that at high densities it is a first 
. ■ order phase transition. In the (J. — T phase diagram, from first chiral phase transition to second order phase transition 
there exists a tri-critical point(TCP). For deconfinement phase transition, it has not a good order parameter except 
for infinite quark mass limit, at which the Polyakov loop severs as an order parameter [13l . [T^l . In recent studies the 
Polyakov l oop has been combined into the chiral models, such as Nambu-Jona-Lasinio modelfla. [l6j and linear sigma 
model |17h19| . which allows to investigate the deconfinement phase transition within the chiral models. Though the 
Polyakov loop is not a good order parameter, it still serves as an indicator of a rapid crossover towards deconfinement. 
As we know in the Landau theory, for the study of the phase transition and the transition order, one should find a good 
order parameter. Once it is identified, the thermodynamic functions could be expanded over this order parameter and 
the transition order could be well studied. For the deconfinement phase transition, besides the Polyakov loop, one can 
also search for other proper order parameters in the effective field models. In the earlier studies of deconfinement, the 
bag models had been often used to investigate the confinement mechanics and the thermodynamics of deconfinement 
phase transition. In this paper we wish to use the effective bag model to study the deconfinement phase transition 
and mainly focus on the study of the transition order and the possible phase diagram of the deconfinement, especially 
the possible influence on the phase diagram by the fluctuations. 

The model we used here is Friedberg-Lee(FL) soliton bag model. The FL model has been widely discussed in past 
decays [20I422T ] . It has been very successful in describing phenomenologically the static properties of hadrons and their 
behaviors at low energy. The model consists of quark fields interacting with a phenomenological scalar field a. The 
a field is introduced to describe the complicated nonperturbative features of QCD vacuum. It naturally gives a color 
confinement mechanism in QCD theory. The model has been also extended to finite temperatures and densities to 
study deconfinement phase transition |23l427| . Here we will try to identify the proper order parameter in this model 
and make an analysis of deconfinement phase transition. 

The organization of this paper is as follows: in section 2 we give a brief introduction of the FL model. The thermo- 
dynamic potential is derived and deconfinement phase transition is discussed at finite temperatures and densities at 
mean field theory (MFT) approximation. In section 3, we make a Landau expansion of the thermodynamic potential. 
In this way the transition order is studied by analyzing the Landau coefficients. By an ansatz of Landau coefficients 
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we discuss the possible phase diagram of deconfinement in FL model. The last section is the summary. 



II. THE THERMODYNAMIC POTENTIAL AND DECONFINEMENT PHASE TRANSITION IN FL 

MODEL AT MFT 



We start from the Lagrangian of the FL model, 



C = jfand" - ga)i> + -(d^a)(d^) - U(a), 



where 
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ip represents the quark field, and a denotes the phenomenological scalar field, a, 6, c, g and B are the constants which 
are generally fitted in with producing the properties of hadrons appropriately at zero temperature. We shift the a 
field as a — > a + a' where a and a' are the vacuum expectation value and the fluctuation of the a field respectively. 
Then the lagrangian becomes 



Ceff = ~ m q )^ + -(VXfl'V) - \rala' 2 - U(a), 
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where 



U(a) 
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m q = ga and vr? a = a + ba + \ca 2 are the effective masses of the quark and a fields respectively. The interactions 
associated with the fluctuation a', such as er' 3 , er' 4 and -tpa't/j, are neglected in MFT approximation. 
According to finite temperature field theory, the partition function is 



Z = 



dip] [da 1 ] exp 



(8 
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where fi is chemical potential of quarks. Completing the integration in partition function Z, together with the 
thermodynamic potential: f2 = —TlnZ, at mean field level, we could obtain 



Q = U(a) + 



1 



d 3 p 

(3 J (2tt) 3 



ln(l 



d 3 p 

(2tt) 3 



ln(l + e - 0(E ->-^) + ln(l + e -^ E «+^) 
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where fi is the inverse of the temperature T and 7 is a degenerate factor that 7 = 2(spin) x 2{ flavor) x 3(co/or). In 
addition, E a — ^Jf?- 



ii and E n = 



p 1 



In our calculation, the parameters are chosen to be a = 17.7 fm~ 2 ,b = — 1457. 4/m -1 , c = 20000, g = 12.16. The 
effective mass of a field is fixed at m a — 550MeV 25j. Then one could plot il versus a for different T as shown in 
FiglU At zero temperature, where f2 = U(a), there are two minima of the thermodynamic potential: one corresponds 
to the perturbative vacuum at a — 0, another corresponds to the physical vacuum at a = a v . The system is stabled 
at the physical vacuum at a — o~ v . It is well known that at this time the quarks are confined in a soliton bag, and the 
system is in a hadronic phase. With temperature increased, the physical vacuum a — a v is lifted up, while the quarks 
has been still confined until the two vacuums degenerate. At this time the deconfinement phase transition occurs, 
and the phase transition temperature is T — T c . After that, the system is stabled at the perturbative vacuum a = 0, 
where the quarks are deconfined and the system is in a deconfined phase. This is a first order phase transition. 

One can also plot the f2 versus a at different n for T = 5QMeV as shown in FigfSJ The deconfinement phase 
transition takes place at [i = [i c where the two vacuums degenerate. The analysis of deconfinement phase transition 
at finite chemical potential is similar to that at finite temperature. 

One can obtain the fi — T phase diagram as shown in Fig|3J In the whole fj, — T phase plane, the transition is first 
order. 
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FIG. 1: The thermodynamical potentials for different temperatures and zero chemical potential: To = OMeV, T\ — lOOMeV , 
T c = 121MeV and T 2 = ISOMeV. 
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FIG. 2: The thermodynamical potentials for different chemical potentials and fixed temperature at T = 50MeV: fii = WOMeV , 
fi 2 = 2QQMeV, ^ c = 255MeV, and ^ 3 = SOOMeV. 



III. A LANDAU EXPANSION AND THE POSSIBLE PHASE DIAGRAM OF DECONFINEMENT 

PHASE TRANSITION 

In above discussion, we know at MFT approximation in FL model the deconfincment phase transition is first order. 
One can plot the a as a function of T, as shown in FigHJ It could be seen that at T = T c , a jumps from nonzero 
value a — a v to zero value <r = 0. In confined phase a ^ 0; in deconfined phase a — 0. Here a could be viewed as an 
order parameter of deconfinement phase transition in FL model, so we can do a Landau expansion of f2 based on a 
and make a thorough investigation of the phase transition order. 

At the MFT approximation, from equation ([5]) the thermodynamic potential could be power expanded by a with 
a 2 , er 3 and a 4 . However, the analytical forms of the coefficients of the expansion are difficult to be obtained. Here we 
will write down the effective form of the expansion as 

n = \A{T, v)a 2 + is(T, ^ + Ic(T, n)a\ (7) 

where A(T, /i), B(T,fj,) and C(T,fi) are the effective parameters which could be determined by a numerical fitting 
process. That means at certain T and /i from the configuration of the f2 versus a one could fit the curve by the 
a 2 , a 3 and a 4 to obtain the values of A(T,n), B(T,/i) and C(T, /x). By the equation (J7J, from Landau theory, it 
is clear that the cubic term a 3 plays crucial role in determination of the transition order. At MFT approximation, 
the fitting results indicate that B(T,fi), as a negative value, will keep decreasing with temperature and/or chemical 
potential increasing. That means this term will never be zero, therefore the transition order of deconfinement at MFT 
approximation can only be first order. 
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FIG. 3: The [i — T phase diagram of deconfinement at MFT in the FL model. 
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FIG. 4: a as a function of T at zero chemical potential in the FL model. 



Now we suppose equation (|7|) is the general form of expansion of thermodynamical potential by order parameter 
a in FL model. And we regard the corrections coming from the fluctuations will effectively modify the parameters 
A(T, fi), B(T,n) and C(T, fj,). In principle they could be calculated by self-consistently resumming the higher order 
loop diagrams led by the fluctuations of a'. However it is very difficult to evaluate these corrections in this way. In 
the following we will treat the coefficients A(T,/j), B(T,/i) and C(T,/x) as the free parameters and make a general 
study of the phase transition order on the FL model by Landau theory. 
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FIG. 5: Phase diagram of deconfinement on the A — B plane in the FL model. 
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In Landau theory, one can make a derivative of the thermodynamic potential to er as 



^=A(T,^ + \b(T,^ 



C(7»ct 3 = 0. 



One can obtain three solutions: 



<7l =0, 



02,3 



-3B ± V9B 2 - 24AC 
2C ' 



(8) 



(9) 



In our case, we assume C > which guarantees that the vacuums are the minima. When 3B 2 < SAC, there is only 
one minimum at a = 0. When 3B 2 > 8 AC, there are two minima. They correspond to the perturbative vacuum at 
(f = and the physical vacuum at a = a v . When the two minima degenerate, one can obtain the condition that: 
B 2 = SAC, at which the deconfinement phase transition takes place. Thus one can draw the critical line of the 
deconfincmcnt phase transition in the plane of B versus A as shown in Fig(5] The phase plane has been divided into 
two parts: the left area beside the line in the plane represents the confined phase, while the right area the deconfincd 
phase. By analyzing the variation of the vacuum, one can obtain that the deconfinement phase transition can be 
either first or second order. If the system goes across the critical line at B ^ 0, the transition is first order. If the 
system goes across the line at B = 0, the transition is second order. 
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FIG. 6: The thermodynamical potentials for different temperatures and zero chemical potential. 
WOMeV, 150AfeV, 180MeV and 2QQMeV from bottom to top. 



The temperatures are 



From above discussion by Landau theory, we know there may be a second order phase transition in FL model, while 
at MFT level, the deconfinement phase transition can only be first order. But if we consider fluctuations beyond 
MFT, there are maybe additional terms which cancel the cubic a 3 term. The second order phase transition may 
be possible. That means the parameter B(T, /i) will go to zero before the transition takes place. The system will 
evolve from left area to right area across the critical line by the axis origin in the FigJSJ In our former calculation 
at MFT, the fluctuations of a' in the Lagrangian have been neglected. These terms are possibly important in the 
cancellation of the cubic term. However, it is very difficult to calculate the thermodynamic potential including these 
fluctuations from the Lagrangian in FL model. In the following, we will make an ansatz based on the form of the 
Landau expansion of the thermodynamic potential to mimic the deconfinement phase transition which have both first 
and second order phase transition. 

We can devise a possible variation pattern of A(T, ^), B{T, fi) and C(T, fi). We suppose at finite temperature and 
zero chemical potential, the absolute value of B{T, n) keeps decreasing and tends to zero with temperature increasing, 
while A(T, /i) first decreases to a negative value and then increases with temperature increasing. C(T, /i) keeps positive 
in all the cases. By this kind of variation, from FigJSJ one could see that the system will evolve from the confined 
phase to the deconfined phase across the axis origin, and the transition will be second order. Thus we make the 
following ansatz of A(T,^l), B(T,/j,) and C(T,n) as 



A(T,n) = a 



(T-T c )(fc 1 T--) + A lM 2 



B{T,fi) = bExp 
C(T,n) = c, 



T + T C 



(10) 

(11) 
(12) 
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FIG. 7: The thermodynamical potentials for different chemical potentials and zero temperature. The chemical potentials are 
350MeV, 392MeV, 420MeV and 470MeV from bottom to top. 

where a, b and c are the parameters of the FL model which have been already given in section II. k\ = 4fm 2 , &2 = 
0.15, Ai = 0.5 fm 2 and A2 = 0.1 fm are the effective parameters of the ansatz. T c is the critical temperature of the 
transition at zero chemical potential which could be seen in later analysis. It also serves as a temperature scaling 
factor which value can be taken as T c = 180MeP. When T = \i = 0, it is clear that A(0,0) = a, B(0,0) = b 
and C(0, 0) = c. One should notice that in our ansatz with the temperature increasing the parameter B(T,[i) will 
be infinitely close to zero but not zero. However when the second order phase transition takes place, the absolute 
value of B(T,ii) will be sufficiently small. At zero chemical potential, from equation (flO|) . one could see at T = T c , 
A(T C , 0) = 0. At the same time B(T C , 0) « 0. Thus the deconfmement phase transition at zero chemical potential and 
finite temperature takes place at T — T c = 180MeP and the transition order is second order. At zero temperature, 
from equation one could see that B(0, fx) will never be zero with chemical potential increasing, which means the 
transition will be first order at zero temperature and finite chemical potential. 

We can also evaluate the thermodynamic potential for different chemical potentials and temperatures. At finite 
temperature and zero chemical potential, the thermodynamic potential as a function of a is plotted in Fig|6] It is 
clear that the phase transition is second order. At zero temperature and finite chemical potential, it could be seen 
from FigJ7]that the transition is first order. The deconfmement phase transition could be presented in a fi — T phase 
diagram as shown in FigJSJ From first order phase transition to second order phase transition there exists a TCP. The 
phase diagram is qualitatively consistent with that of the chiral phase transition. However, how to obtain the credible 
phase diagram of deconfmement through the direct calculations including the fluctuations from the Lagrangian of the 
FL model deserves a further investigation. 
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FIG. 8: fi — T phase diagram of deconfmement with TCP in the FL model. 
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IV. SUMMARY 

In this paper we have discussed the possible phase diagram of deconfmement in FL model. By the calculation only 
in the MFT approximation and without the fluctuations, the deconfmement phase transition can only be first order at 
finite temperature and chemical potential. By the Landau expansion of the thermodynamic potential and the analysis 
through Landau theory, we show that the deconfmement phase transition can also be second order, which will not 
appear in the MFT approximation but will possibly appear when nonlinear fluctuations are considered. Thinking of 
the difficulties in calculating the fluctuations, we have not done the calculation here but made the ansatz that the 
Landau coefficients are certain functions of temperature and chemical potential. By this ansatz we obtain the possible 
fJ. — T phase diagram of deconfmement in FL model which is similar to that of the chiral phase transition. That means 
the deconfmement phase transition is first order at low temperature and high chemical potential while second order 
at high temperature and low chemical potential. From first order to second order phase transition there exists a TCP. 
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